Sharp interface limits of phase-field models 
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The use of continuum phase-field models to describe the motion of well-defined interfaces is discussed 
for a class of phenomena, that includes order/disorder transitions, spinodal decomposition and 
Ostwald ripening, dendritic growth, and the solidification of eutectic alloys. The projection operator 
method is used to extract the "sharp interface limit" from phase field models which have interfaces 
that are diffuse on a length scale ^. In particular, phase-field equations are mapped onto sharp 
interface equations in the limits <C 1 and £,v/D <^ 1, where k and v are respectively the interface 
curvature and velocity and D is the diffusion constant in the bulk. The calculations provide one 
general set of sharp interface equations that incorporate the Gibbs-Thomson condition, the AUen- 
Cahn equation and the Kardar-Parisi-Zhang equation. 
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I. INTRODUCTION 

Many inhomogeneous systems involve domains of well- 
defined phases separated by thin interfaces. These in- 
clude non-equilibrium systems where phases are separat- 
ing by spinodal decomposition or nucleation and growth 
[01 , or where solidification occurs by dendritic growth |^] 
or by the growth of eutectic crystals j|] . The phenomeno- 
logical description of these phenomena involves the mo- 
tion of well-defined sharp interfaces. The origin of such 
descriptions is often transparent, being obtained by sym- 
metry arguments and common sense. Nevertheless the 
properties of sharp interface models can be quite subtle 
as is the case for dendritic growth. 

Unfortunately sharp interface models are difficult to 
simulate since this usually involves solving a diffusion 
equation subject to moving boundary conditions at the 
interfaces. A more convenient approach is to simulate 
models which describe the bulk phases as well as the 
interface structure. While these models are wasteful in 
terms of simulating bulk regions [Q , no explicit boundary 
tracking is needed. This is the key element to a popular 
method for studying systems out of equilibrium, called 
"phase-field" modeling. In such an approach one or more 
continuous fields which are functions of space r and time 
t, are introduced to describe the phases present. Typi- 
cally these fields vary slowly in bulk regions and rapidly, 
on length scales of the order of the correlation length, ^, 
near interfaces. The free energy functional determines 
the phase behavior and, with the equations of motion, 
gives a complete description. In other contexts, such as 
critical dynamics 0J^, the fields are the order param- 
eters distinguishing the phases. In a binary alloy, for 
example, the local concentration or sublattice concentra- 
tion can be described by such fields. The ideas involved 
in this approach have a long history, going back to van 



der Waals [|6|. Within the materials community, the use 
of continuum field models is associated particularly with 
the work of Cahn and collaborators 0,^. 

Phase field models provide a complete microscopic de- 
scription, but are not necessarily appropriate for a par- 
ticular system. These models apply to a large number of 
microscopic systems only when ^ is much larger than any 
particular microscopic length, such as the lattice spacing, 
involved in the surface structure. Interpreted in the sense 
of describing the universal features of many microscopic 
models, f is a mesoscopic length representative of mi- 
croscopic structure. In a similar manner sharp-interface 
models apply to a large number of continuum field mod- 
els of pattern formation or phase separation in the limit 
that the length scales defined by the patterns are much 
larger than ^. An important difference exists however in 
the construction of the two approaches. Standard sharp- 
interface models are constructed from phenomenological 
descriptions of interfaces, while phase field models can be 
constructed to explicitly obey the fundamental principles 
of statistical mechanics. 

While continuum phase field models provide a funda- 
mental approach which is clear and workable, it is impor- 
tant to establish the connection between this description 
and the sharp interface description. The main difficulty 
which arises is how to take account of the finite thick- 
ness ^ of the diffuse interface of the continuum model. 
There has been a great deal of discussion in the literature 
on how this process is to be undertaken. Some workers 
have extracted the interface equations by taking the limit 
where the interface width of the phase field model goes to 
zero Q for the Stefan problem of a pure material. This 
approach is not very useful since the interface width is 
always finite. More recently these calculations have been 
extended for special choices of the free energy functional 
J- to include an interface of non-zero width 
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The purpose of this paper is to provide a clear and 
detailed derivation of the sharp interface equations. 
The approach follows the projection operator method 
of Kawasaki and Ohta is generally applicable and 
eliminates the counter intuitive necessity of limiting the 
derivation of the sharp interface limit to some artificial 
form of the free energy functional T [|lO|. This general 
calculation provides one set of equations that relate the 
parameters of the phase field equations to those of the 
sharp interface equations for a broad class of phenomena 
including order-disorder transitions, dendritic growth, 
phase separation in binary alloys, eutectic growth, and 
surface roughening. In particular, thermodynamic con- 
sistency is automatic in the present approach for nonzero 
interface widths. This is in contrast to the unphysical 
approach of taking the limit of a zero interface width, 
which requires fine tuning of the free energy to obtain a 
thermodynamically consistent theory. 

The structure of the paper is as follows. In the next 
section a general introduction to phase field and sharp in- 
terface equations is given. This is followed by a detailed 
calculation of the sharp interface equations from a set of 
phase field equations. The usefulness of the sharp inter- 
face equations is then illustrated by considering the linear 
stability of planar interface separating various phases. 



II. DESCRIBING INTERFACE PHENOMENA 
USING PHASE-FIELD MODELS 

To sec the relationship between the two approaches, it 
is useful to construct the simplest equilibrium and non- 
equilibrium descriptions of surfaces by both methods. 
The main ideas developed below appear in the more gen- 
eral case discussed in the following sections. These ideas 
are present in work by Allen and Cahn [|| on the mo- 
tion of anti-phase boundaries in the kinetics of an order- 
disorder transition in a binary alloy. This work recov- 
ers the sharp-interface description from the phase-field 
model in the appropriate limit and predicts an experi- 
mentally testable consequence of a finite ^. 

First, consider a surface in equilibrium. The reason 
why droplets are spheres and interfaces are locally flat is 
that these shapes minimize surface area. The excess free 
energy of a surface is proportional to its area A 



l^T = aA 



(1) 



where the proportionality constant u is the surface ten- 
sion. This simple common sense approach is the sharp 
interface model. 

Contrast this with the phase field approach. The free 
energy functional J-'{'ip) for the scalar order parameter 
ipi'f^ has the following familiar form, consistent with re- 
flection symmetry ^-{'ip) = JF(— -0) 



dr 



rl 



(2) 



where > so that the square gradient interaction 
gives the free energy cost for inhomogeneities and the 
local bulk free energy has a double well form satisfying 

f{ip) = f{—ip). An example is the form 



/ 



(3) 



where a oc (To — T) > with To the mean field critical 
temperature and 6 > is a constant. 

The motivation for employing such free energy func- 
tionals is as follows. First, the free energy functional 
is constructed of local or bulk terms, f{ip), which in- 
teract through the gradient term. Secondly must be 
an analytic function since it describes a microscopic sys- 
tem and should not be confused with the true thermo- 
dynamic free energy, F = — fcsT In^^ exp(— .?^/A:bT) 
where /cb is Boltzmann's constant, which is non analytic 
at a phase transition fl^ . The square gradient term is 
the simplest way for the model to be well defined on 
small length scales. Higher-order gradient terms, such 
as K4(V'^tp)'^ could, in principle, describe correlations 
on still smaller scales, such as the internal structure of 
the interface. This is analogous to what is described by 
the Percus-Yevick or hypernetted-chain theories of dense 
fiuids |T^. In the absence of any specific experimental 
motivation, however, terms proportional to K4 may be 
ignored. From the point of view of identifying and study- 
ing a well defined microscopic model, the square gradient 
free energy provides a complete description. Alterna- 
tively, the square gradient theory describes a large class 
of microscopic models provided we consider mesoscopic 
scales £,{K)/ao ^ 1, where Qq is a lattice constant. 

It is particularly convenient that the microscopic phase 
field description is so "close" to mean field theory. In 
that approximation, / is simply the bulk free energy den- 
sity. In practice, the form of the free energy functional 

can be constructed straightforwardly with reference to 
the phase diagram of interest. While one can also con- 
struct microscopic lattice gas models of phenomena such 
as phase separation and dendritic solidification ||l5|] , us- 
ing similar arguments of universality and simplicity, such 
models do not have this convenient feature. 

If fluctuations are small, the equilibrium behavior of 
the model is determined by the mean field approximation 
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dl 

9V 



0. 



(4) 



The homogeneous bulk solutions, valid well below Tc, are 
given by df/dipeq = and are equal to. 



-ip = ±'4>eq = ±\/a/b. 



(5) 



for the bulk free energy given in Eq. (|3|). Fluctuations 
around the bulk solutions satisfy 



(6) 
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for r ^ ^, where ^ = K,p/2a and (• • •) represents an 
ensemble average. 

Now consider a system with a flat interface located at 
y — Q when ^{r) depends only on y and the interface 
profile il)™{y) is the solution of 

with ^'"(±00) = dzipeq- Solving by quadratures gives 

(8) 

from which Eq. (|^) follows immediately. In contrast to 
the sharp interface approach, this yields an explicit form 
for the surface tension 



(Tw, = dy 



/dip" 



\ dy 



(9) 



For the particular form of ^{1^) of Eq. (g) and Eq. (|^), 
one finds 



V''"(y)=^e,tanh(^ 



and 



3C 



(10) 



(11) 



The contrast between the sharp-interface model and the 
phase- field model is revealing. Both give the same macro- 
scopic description but the phase field approach gives a 
complete description down to the smallest length scale. 

Now consider a sharp interface that is in local but not 
global equilibrium due to a gentle curvature. For sim- 
plicity consider phenomena where the order parameter is 
not conserved, such as occurs when a binary alloy under- 
goes an order-disorder transition and ip is the sublattice 
concentration. The interface moves to locally reduce the 
surface area and surface free energy with an interface ve- 
locity V. Expanding u in a Taylor series in powers of k 
gives 

V = -VK (12) 

to lowest order in the curvature. This is the sharp in- 
terface theory for the motion of anti-phase boundaries. 
Note that, since v is the only coefficient which enters the 
theory and has dimensions of a diffusion constant, any 
time-dependent length R{t) must satisfy 



R{t) = [utr, 



(13) 



by dimensional analysis, where the growth exponent 

n=l/2. (14) 



Such an approach was first done by Lifshitz | [T^ and 
by TurnbuU ||l^. The sharp interface treatment alone 
cannot predict the value of v. An additional argument, 
which turns out to be incorrect for the motion of an- 
tiphase boundaries, was used to predict u oc a^T^^, where 
is a mobility. 

A first principles approach to this phenomenon is due 
to Allen and Cahn g]. Neglecting noise, the equation of 
motion for the non-conserved sublattice concentration is 

Sli 



dt ~ *6ip ~ ^ 



(15) 



Allen and Cahn denoted the position of the antiphase 
boundary by a curved, time-dependent interface M(r, t) = 
0. They then looked for solutions of the form ip{r, t) = 
V''"(w(r,i)). This gives 



dip _ c'V''" 
dt du 



dv? 



K^K— h — 

ou 



a/1 

dip. 
(16) 



where k — — V • n with n ~ Vw/|Vu| the unit vector 
normal to the interface and u{f,t) a coordinate in the 
direction n. Eliminating df /dip™ using Eq. (|^) gives 
Eq. (12) where, in contrast to the sharp interface the- 



ory, one obtains an explicit expression for the transport 
coefficient. 



(17) 



Results in the presence of stochastic noise have been 
obtained by many authors, particularly Bausch, Domb, 
Janssen and Zia ]l^ , and by Kawasaki and Ohta ||l^ . 

Although both approaches correctly find that the nor- 
mal velocity is proportional to the curvature, the AUen- 
Cahn result for v is noteworthy. The earlier theory, 
which argued v cx a^T^ implies a strong dependence 
of the velocity v on interfacial thickness since a cx 1/^, 
from Eq. (pi]). In contrast, Allen and Cahn predict that 
v is independent of interfacial thickness. This was clearly 
demonstrated in an experiment by Pindak et al. |l9| ] 
where they studied orientational patterns in freely sus- 
pended dipolar smectic C liquid crystal films. Since their 
smectic C films have a permanent electric dipole moment 
of magnitude P, the director angle (p can be oriented with 
an electric field of magnitude E. The free energy is 



df 



-K{V(pf - EP cos^ 



so the width of the interface, 



^ cx 1/VS, 



(18) 



(19) 



can be varied easily. The experiments directly verify 
that the size of a domain of stable orientation grows as 
R{t) = (^^^)l/^ where V is independent of interface width 
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in accord with the prediction of Allen and Cahn Q . In 
addition, the experiments show that the sharp interface 
result, TJ oc K, is independent of interface width provided 
1. 

It is straightforward to include an external field c cou- 
pling linearly to ip. Then T ^ T + Text with 



(20) 



and one of the phases becomes metastable depending on 
the sign of c, so that the interface translates even if it is 
flat. Hence 



(21) 



with vq (c) cx c. This allows one to simulate Kolmogorov- 
Avrami-Johnson-Mehl growth of droplets |2^] and, in the 
presence of noise, Kardar-Parisi- Zhang dynamic rough- 
ening ||. 

Other processes can be simulated when one of the 
phases is metastable and the growth of the stable from 
the metastable phase is controlled by a conservation law. 
This describes solidification of a metastable supercooled 
liquid phase and the growth of the stable solid is limited 
by the diffusion of latent heat from the surface of the 
moving solid front. The external field c is then propor- 
tional to the latent heat. In the sharp interface formula- 
tion, c obeys a diffusion equation in the bulk phases 



dc 
dt 



(22) 



where Dc is a diffusion constant. 

The steady state velocity v of the interface is given by 
integrating Eq. (p2|) across the interface to obtain 



{DtVc\+ - D-Vcn 



(23) 



where the superscripts ± refer to the values of the nor- 
mal gradient of c on either side of the interface. The 
condition of local equilibrium at the interface is 



Sc 



(24) 



which is a Gibbs-Thompson condition relating the local 
excess concentration Sc to the curvature. This says that 
excess external driving force is balanced by the curva- 
ture K. In regimes of high undercooling, this is some- 
times supplemented by an additive term proportional to 
V, describing kinetic undercooling. 

To study this by a phase-field approach, let c be a 
continuous function of space and time which is conserved: 



dc (r, t) 
dt 



n5T 



(25) 



where Fc is the mobility of the field c. For the model to 
be well defined, a self energy for c or a positive additive 
contribution to T of the form 



j dr 









(26) 



must be included. Within a mean-field approximation, 
this gives a homogeneous equilibrium solution c cx ?/;. 
Note that the interface invades the metastable phase be- 
cause of Text but it must also satisfy the conservation of 
c as defined by Eqs. (25) and (p6|). This implies that the 
interface deforms into a parabolic shape, dumping excess 
c to the sides while propagating forward at a constant 
velocity. The parabolic shape has a constant growth ve- 
locity in the forward direction, satisfying the fact the 
system is driven with a constant thermodynamic force, 
while lateral growth has a velocity ~ t^^/^, thereby sat- 
isfying the conserved diffusion equation for c . 

It turns out that, although this is the right approach, 
the implementation needs some fine tuning. First, when 
considering dendritic growth, the theory of microscopic 
solvability ||2^] has shown that dendrites require an 
anisotropic surface tension to be well defined. Hence, 
one must let K — > K (\J if) / \V ^\) in some convenient pre- 
scribed way. Next, in a very useful paper, Kobayashi 
p3| has noted that, to keep the equilibrium solutions ij^eq 
from shifting appreciably when c is applied, Text of Eq. 
(^ ) should be modified to 

Text = - J drc-^iiP) (27) 

where 'I'(V') is an odd function of ip satisfying 

d^/dyJeq = 0. (28) 

For example, if i/jeq — ±1, one can choose d'^/dtp — 
(1 — V'^)^- Other forms are possible. 

Finally, one can choose to make the c field phase sep- 
arate by replacing the self-energy term in the free energy 
of Eq. ( |26| ) with a double well form, analogous to Eq. (|^). 
This permits the study of eutectic crystallization In 
the remainder of this paper a detailed implementation 
of these ideas is given, making connections to the sharp 
interface limit. 



III. EQUILIBRIUM 

Consider two fields, a non-conserved phase field ip and 
a conserved field c. The phase field distinguishes be- 
tween, for example, liquid and solid phases and the c 
field can be taken as a concentration. The free energy 
functional describing the system can be written as 



T{^,c} 



dr 



(29) 



where /(■)/), c) is the local bulk free energy density and 
the gradient terms account for interfaces and inhomo- 
geneities as discussed above. The dynamics of these fields 
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are described by the equations of motion for the non- 
conserved tp, Eq. and the conserved concentration 



c, 



oc 



(30) 



where 5J-{^l),c\/5c — ^ is a, chemical potential and Tc is 
a mobility. The usual additive noise terms, related to the 
transport coefhcients by fluctuation-dissipation re- 
lations , have been neglected for simplicity. In mean 
field theory, the equilibrium states of the system with an 
interface normal to ij are defined by 



5/ 



and 



dl 

dc 



fJ-e 



(31) 



(32) 



where /ie<j is the chemical potential of the uni- 
form equilibrium states. Integrating Eq. (^) over 
c gives the Maxwell's equal area construction rule, 
J^^ dc {df/dc~iieq) — 0, where ci and C2 are defined 

by 9//(9c|ci,c2 = Meg- 

For a generic /(c, ip) there may be many possible equi- 
librium and metastable states contained in this free en- 
ergy. For illustration, consider the following bulk free 
energy 

/(V', c) = [c In c + (1 - c) In (1 - c)] 

J A/ 

+ (^aAT - /3(c - i)^) ^{yj) - + (33) 

where *(V') = 2^^ - 4V'V3 + aV'VS, AT = {T-Tm)/Tm 
with Tm the melting temperature and the other phe- 
nomenological parameters are determined by matching 
to experimental phase diagrams. If these parameters are 
chosen a.s a = /3 = 1.0 and u — 0.6, the mean-field phase 
diagram shown in Fig. (|l|) emerges. As can be seen, 
this phase diagram contains liquid/solid and solid/solid 
coexistence regimes. For this symmetric free energy the 
melting temperature at c = 1/2 is denoted Tm and the 
critical point of the solid/solid coexistence regime is at 
(c,r) = (1/2, Tc) with Tc < Tm- As the parameter 
u is decreased, Tc increases until the solid/solid coex- 
istence region collides with the liquid/solid coexistence 
regime when Tc > Tm and a eutectic point is formed at 
(c,T) = (1/2, Tb) as shown in Fig. (|). 

As can be seen from the phase diagrams of Figs. (|l|) 
and (|2|), this simple free energy contains many phases 
and, in conjunction with appropriate equations of mo- 
tion, can be used to study a wide variety of phenomena. 
A number of 'quenches' have been highlighted on these 
diagrams to illustrate several different kinetic processes 
that may arise. 
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FIG. 1. Mean field phase diagram obtained from the bulk 
free energy of Eq. (M) for the parameters a = (3 = 1.0 and 
u — 0.6. In this figure the regions containing vertical and 
horizontal lines are liquid/solid and solid/solid coexistence 
regions respectively. 




FIG. 2. Mean field phase diagram obtained from bulk free 
energy given in Eq. (^3|) for the parameters a — [3 — 1.0 and 
u = 0.45. In this figure the regions containing vertical and 
horizontal lines are liquid/solid and solid/solid coexistence 
regions respectively. 

A quench is defined as a rapid change in temperature 
which takes a system from one region of the phase dia- 
gram to another and is often considered instantaneous in 
theoretical modeling. In the next section the dynamics of 
an interface separating a stable and a metastable phase is 
considered. The calculations are done in a general man- 
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ner to include all the possible quenches shown in Figs. 
(H) and (H). To construct the relevant sharp interface 
equations no reference will be made to the explicit form 
of the bulk free energy term /. It is assumed that / has 
been chosen merely so that all the phases of interest are 
well defined. 



length scale ( such that l<<C/^<<l/e, then distinct 
'inner' and 'outer' regions can be defined, as shown in 
Fig. and it is appropriate to solve in both inner and 
outer regions and match the solutions at the length scale 
Formally, the technique requires an inner expansion 
near the interface and an outer expansion far from the 
interface. 



IV. NON-EQUILIBRIUM 

The goal of this section is to derive the sharp interface 
equations for a systems described by a free energy func- 
tional (such as is given by Eqs. (|29|) and ( p3[)) and the 
Langevin equations given in Eqs. ([15[) and (|30|). These 
latter equations should be supplemented by additive 
stochastic noises, 7y^(r, t) and r]c{r,t) of zero mean and 
correlations < r]^{r,t)r]^{f',t') >= 2V^T5{r-r')5{t-t') 
and < ric{r,t)r]c{r\t') >= 2TcT\/^S{r^ f')S{t - 1') with 
< VipVc >= as required by the fluctuation dissipation 
theorem. With these stochastic noises, the dynamical 
equations Eq. (|l5| ) and Eq. (|3^) are the simplest equa- 
tions which respect the macroscopic conservation laws 
and also ensure that the system evolves towards its ulti- 
mate state of thermal and mechanical equilibrium with 
its external environment. 

The asymptotic analysis proceeds by expanding 
around a planar equilibrium interface described by Eqs. 
(|l|) and (||). The interface can be taken out of equilib- 
rium by either gently curving it or by making one of the 
two bulk phases metastable. In the former case, a gen- 
tle curvature is one in which the radius of curvature l/n 
is large compared to the interface width or correlation 
length. Thus one small dimensionless expansion param- 
eter is K^. In the latter case, the difference between the 
free energy of the stable and metastable phases causes 
the interface to propagate into the metastable phase. If 
the free energy difference is small the propagation veloc- 
ity V is small. In this context, a small velocity means that 
the interface moves so slowly that a steady state diffu- 
sion field is allowed to form in front of the interface. In 
other words, the time for the diffusion field to relax when 
the interface moves a distance ^ should be much smaller 
than the time £^/v taken for the interface to move that 
distance. Since the diffusion time r — ^'^/D, this leads 
to another small dimensionless parameter ^v/D which 
is known as the interface Peclet number. In the follow- 
ing analysis the interface equations will be obtained to 
lowest order in both small parameters. Technically the 
expansion to lowest order in both small parameters can 
be achieved if they are regarded as the same order in the 
expansion. In the calculations to follow both parameters 
will be taken to be 0(e) with e << 1. 

The calculations make use of the fact that the fields 
behave qualitatively different close and far from the in- 
terface. In the region close to the interface, the fields 
vary rapidly over distances 0{^) while, far from the in- 
terface, they vary over distances ©(^/e). If there exists a 




c(u) 



FIG. 3. Illustration of inner and outer regions used in com- 
putations. 



A. Inner Expansion 

Consider an inner region defined by —C. < u < +(^, 
where u is a coordinate normal to the interface and 1 < < 
C/^ << 1/e. The aim is to obtain asymptotic expansions 
for the solutions to the evolution equations (|l^ and ( |30| ) 
valid in this inner region. The latter can be written in a 
compact form as 



1 <9c „2r 



where 



df 

6fl = ^i{f) - lleq = -Kc\/'^c + fl 

oc 



(34) 



(35) 



The first step is to partition the system into two regions 
V-^ and bounded by surfaces S+ and S- respectively. 
The region is defined by < it < ^ and and similarly 
for The position of the interface between two bulk 
regions is defined as u{f,t) = 0. These definitions are 
purely formal, but to fix ideas, the surface u{r, t) = may 
be regarded as the surface near which the fields c, ^ vary 
rapidly over distances C'(^). It is then useful to define 
Green's functions G^{f,f') in the regions — C < w < 
(G") and in < u < C (G+) obeying 
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(36) 



and satisfying the boundary conditions, G{r,r') — 
at w = and u' — 0, dG{r,r')/du = at u = 
dG{r,'r')/du' = at u' — ±C, and periodic in the other 
directions. Note that both r and r ' he in the same region, 
V+ or V-. 

Multiplying Eq. (^) by G"^ and integrating over f ' e 
V± gives 



dr 



., G^{f,r') dc' 
'dt 



JV± 

j dS' ■ (^Sfi'V'G^ - G^V'S^i' 



(37) 



where Sfi' — Sfi{r") defined as in Eq. ([35[), Jy dr' de- 
notes integration over V+ defined by < u(r ') < C, and 
is+ denotes integration over the boundaries S+ en- 
closing V+ and similarly for V- enclosed by S-. Multi- 
plying Eq. (|3^) by dc^du where c™(m) is determined by 
the solution of Eqs. (|3l|) and (^) for a planar interface in 
thermal equilibrium and integrating over —C, < u < +(^ 
gives 



r+C aJn 
B + S= du^Sfi 



C 



du 



+^ fdf 



du^ ( ^ - JsTrV^ 



_^ OU \ OC 

where B = B+ + B- and S = S+ + 5" with 



B^ 



±- 



1 



c Jo 



du 



du 



dr'G^{r,r') 



v± 



dt 



(38) 



(39) 



± 



du^ j> dS' ■ ( S^i'VG^ - G^^'Sfi' 



s± 



(40) 



An analogous formula for tp is obtained by multiplying 
E£^ (111) by with Vi'H") the solution of Eqs. 

( pi| ) and (|3^), and integrating over —( < u < to 
obtain 



1 



_f du dt 



, d%^ ( 2 df 



(41) 

Each term in the above equations can be systematically 
expanded in powers of e. In this paper, attention is re- 
stricted to the terms 0{e) as much of the relevant physics 
is apparent at this order. Going to higher order in e does 
not yield any new physical insight but does require con- 
siderable book-keeping skill. 

To facilitate the expansion, c{r,t), '4'{r,t) and the 
chemical potential /x(r, t) are expanded in a power series 
in e. 



c(r-,t)=4Xr))- 



+ e'Sd^ + ■ 



e5c 



(42) 



where the superscript 'in' refers to the inner solution. 
To expand the Laplacian in powers of e, it is useful to 
introduce a curvilinear coordinate system with one co- 
ordinate u along the local normal to the interface and 
(d — 1) coordinates s perpendicular to u and tangent to 
the interface. For simplicity a two-dimensional system is 
considered where s is the scalar arc length. Note that the 
©(e") terms, Cq" and V'o" in Eq. ( ^2|), are the equilibrium 
planar interface solutions of Eq. (pi]) and Eq. (|3^) . 

At this early stage of the calculations it is worth point- 
ing out that the exact position of the interface has not 
been specified. The choice of the exact interface position 
is somewhat flexible to within a distance ^ and there is 
no particular reason for choosing the interface position to 
be defined by ^{u{r,t) = 0) = as is often done in the 
literature. Indeed, this particular choice can lead to un- 
reasonable constraints on the free energy for a mapping 
between the phase field model and the sharp interface 
limit to be possible. In this work the interface will be 
chosen to be the Gibbs surface defined so that the ex- 
cess surface concentration is equal on both sides of the 
interface. This is in essence a solvability condition which 
ensures that the chemical potential is continuous across 
the interface. 

The transformation from Cartesian to curvilinear co- 
ordinates (see Appendix A) leads to the formal expansion 



(43) 



where the specific form of £„ depends on the expansion. 
In the inner region, derivatives of the fields with respect 
to u are much larger than derivatives with respect to 
s which are identically zero when the curvature k and 
the Peclet number vanish. This is taken into account by 
introducing the dimensionless variables u and s which are 
0{e^) hy u — S^u and s — ^s/e. As shown in Appendix 
A, this scaling leads to 

^0 = duu 
Ci = fidu 
^2 = dss - K^udu 

where the dimensionless curvature, k = ^K/e, is of order 
unity. 

Lastly, the time derivatives in Eqs. (^) and ( p| ) must 
be expanded in e. For these calculations, it is convenient 
to work in the frame in which the interface is stationary 
so that 



d_ 

dt 



d_ 

dt 



U • V 



(44) 



where v is the interface velocity which has components 
normal and tangential to the interface. The time deriva- 
tive on the right hand side corresponds to relaxational 
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dynamics not accounted for by motion of the interface. 
When this operator acts on the fields, c and ip, the tan- 
gential component and time derivative are of order e"^ and 
can be dropped Thus to ^ 0{e), d/dt\p becomes; 



di 



T OU 



(45) 



where the normal velocity has been expanded in a power 
series in e, 



^ _9u ^ m 

^normal — 777 — / , ^ ^r. 
Ot T 

m— 1 



(46) 



Using these expansions and expanding / around Cq" 
and V'™ the right hand sides of Eqs. ( ^ and (^) become 



du—— 

_^ OU 

_f OU 



(1,0) 



(47) 



and 



f du- 

_^ OU 



+ 



(48) 



where /f"^") = a"+'"//9c"5V'"|^in,,in, if^ = i^c/C', 
i^^ = K^/S."^ and C = C/C Terms" of 0(6°) vanish by 



construction. For later use, it is convenient to perform 
)m 



to pc 

partial integrations on combinations of terms ||27|] 



(2:0) 



du 



-C 



duScTfl 



(1,1) 



du 



(49) 



and 



du 



i9u 



9m ' 



- P 



-C 



du 



(50) 



since derivatives of c™ and -0™ vanish at u = ±C in the 
limit C = C/C » 1- 

To complete the calculation, the left-hand sides of Eqs. 
(^) and (^Tj) are expanded to lowest order in e. The 
expansion for ip in Eq. (^) is straightforward 



^ ' a^ 



where 



fflU „_ t'l- 



9w du 



K 



'du('-f 
-C V du 



(51) 



(52) 



The equivalent expansion for c is more complicated. 
The algebra is given in Appendix B. Formally, the results 
of these calculations can be written as 



B + S = e{Bi/T + Si) + e2(62/T + ^2) 



(53) 



where B„ and 5„ are given in Appendix B. 

Putting all these results together gives the following 
two equations to ©(e) 



ViS, a. 



(54) 



and 



Ac<5Mr(0,s) = -^+Ai 

f^'du[c°,-\u)-ci^{u)r 



C 



du 
dS^i^' 



du 



-CJ-c 



du[cr(«)-ci,"(u)] 



+CJo 



(55) 



duK^\u)-c^iu)] 
where Ac = c™{C) — Co"(— C) is the miscibility gap, 

cr(-)-(jSJ (56) 



K, 



'du('-^ 
du 



-C 



(57) 



and 



A,,r^^(^^^?$.s^?g]!r (58) 



du 



du 



Equation (55) gives the chemical potential fi of the in- 
ner solution at the interface (i.e., at u = 0) which must 
be matched to the outer solution at u = ±C- An expres- 
sion for Sfi^i{±Q can be obtained by expanding Eq. ( |3^ ) 
to 0{e) 
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(59) 



Integrating Eq. ( |59| ) twice, first from w to C and then 
from to C yields 



du 



2 /-C 



(60) 



and similarly for 6^™{—C,). From Eqs. ( p5| ) and ( pQ ) we 
obtain 



Ac(5Aii"(±C, s) = -cTcK/e + Ai ± AcC 



du 



Vie 



2 



rfz2[cr(u)-ci,"(i2)]^ 



du 







d^2[cr(^i)-4"(^)] 

^ ds[cr(«)-c;,"(u)] 

+Ac^ I du[cS"*(fi)-4"(w)] 



du 



(61) 



c Jo 



The integrals in Eq. ( |6l| ) can be written in a more useful 
form by noting that C ^ 1 in the inner region so that 
d^{±u) = 4"(±oo) for |m| > C- Eq. iM) becomes 



l'^^duK{u)-cnu)r 



du 
ddijif 



— OO 





du[cl^\u) - d^(u)] 



du 



+ 00 



du[cZ-\u) - ci,"(u)] 



-Ac^^ / du[cr{u)-d^{u)] 



(62) 



One last result will be needed and is obtained by inte- 
grating Eq. (pg) over —C, < u < 



Ac^^K du 



dSnf 



+C 



du 



-C 



(63) 



The solution for Sii™{±() must be matched to the solu- 
tion in the outer region. 



B. Outer Expansion 

Far from the interface, the fields and c vary slowly in 
space and are close to the bulk equilibrium values ipeq and 
Ceq. Variations of the fields in the bulk regions far from 



the interface take place on length scales 0{^/e) in all di- 
rections which implies that suitable dimensionless space 
and time coordinates are {u,s,t) = (eu/^, es/^, e^i/r). 

Expanding ipif) about the bulk equilibrium solution 
^p{r) — (5-0°"* (f^) -I- 2peq gives 



dt 

1 d^f 



dl 

dip 



5ip° 



2! di}}^ 



1 av 



eq 



eq 



(64) 



By definition, {df /djjAeq = and, since (5-0°"* = in the 
limit e — > 0, Eq. (p4) is linear at 0{e). Furthermore, 
[d'^ f / dip"^) eq > so that ^0g"* vanishes exponentially 
with time for all wavelengths. Thus, (50q"* is trivial in 
the outer region and can be ignored. It is convenient to 
expand c°"* and in the outer region as 



out /^\ ^out 



°"'(f) 



(65) 



where eg"* is given by Eq. ([s^). At 0{e^) in dimension- 
less variables, in the lab frame 



di 



(66) 



where V = (^/e)V is the scaled dimensionless derivative 
suitable for the outer region. This simplifies to a linear 
diffusion equation for the chemical potential inside the 
bulk phases which reads, in dimensional units, 



d^j° 



dt 



(67) 



where = Tc{dfi°^^ /dc°'^^)eq is a diffusion constant. 
The value of Dc depends on the bulk equilibrium phase 
considered. 



C. Matching and the Gibbs Surface 



To solve the diffusion problem of Eq. (|6^), initial val- 
ues 5/if"*({2 — 0, s) are required. These are to be obtained 
by matching to the inner solution 5/x'"(u, s) at u — ^. 
To obtain (5/x?"*(w = 0±,s) from 5/x?"*(w = ±C, s) with 
C = eC, it is useful to Taylor expand about u = ( 



s^l'^'iu,s)^s^,r\±C,s) + {uTC) 



(68) 



In the outer region, C ^ 1 and this expansion is valid at 
u = 



du 



±C 



(69) 



Since 5jj°^\±0 = <Smi"(±C), we can use Eq. (g) and 
Eq. ( |69| ) to obtain 



9 



du 



-CJ- 



ds[cr(«)-4"(u)] 

+CJq 



du 



(70) 



which gives the appropriate boundary value of 
6fi°^^{0,s). The inner solution (5/j,']"(0) differs from the 
outer solution since the matching is done at 



u — and extrapolated linearly to u = by Eq. ( |69[ ) 
The extrapolation and matching of to S^°^^ is i 
trated pictorially in Fig. (Q). 



lus- 



c-(~) 





du |_cg"Hu)-c|,"(u)] 



/ ° du Lci,"(u)-o-'(u)] 




(u) (dashed line) with (5/i^'(^t) 



FIG. 4. Matching of Sfif 
(solid line) at u = ±(^. 



The right hand side of Eq. ([70|) appears to depend 
on whether the inner and outer solutions are matched 
at M = +C or u = —( because of the last term. This 
ambiguity is eliminated by defining the interface to be a 
Gibbs surface at w = defined by the condition 



+ C30 



du[cg"*(M) -cj)"(u)] = 



(71) 



This can always be satisfied by choosing the position of 
the interface at w = to be such that Eq. ([ti] ) is satisfied. 
In essence, the interface position is determined by the 
condition that the excess concentration on one side of 
the interface is exactly compensated by the deficit on the 
other, as shown in Fig. (^. This can be regarded as a 
solvability condition. 



FIG. 5. The Gibbs surface at u = defined by Eq. (|7l|) 
which matches the concentration deficit on one side with the 
concentration excess on the other. 

Another result that will be needed is obtained by 
matching the first derivative of the chemical potential 
is obtained from Eq. (|6^) which reads, in dimensional 
units 



du 



(72) 



Matching derivatives of the inner and outer solutions for 
/i and extrapolating back to m = 0* by Eq. (68), gives 
the standard result 



^'( 



du 



0+ 



du 



(73) 



since /i°^*(it) is linear for < |?£| < ( from Eq. (|6g 
Finally, combining Eqs. ([70|), ( [7l| ) and (|7|), gives the 
chemical potential at a moving, curved interface 



Ac(Ai™*(0,s)-/Xe,) =-(TcK 

where Ai is given by Eq. (|5q) 



+ £^v + Ai+Oie^) 
and 



(74) 



1 



'du[[cr{u)r~K{u)r). (75) 



V. SUMMARY OF RESULTS 

All the results can be combined into a single set of 
boundary layer equations that can describe many differ- 
ent physical phenomena. Typically the boundary layer 
equations are written in terms of the concentration, 
which is the outer region is simply related (to order e) to 
the chemical potential by the relationship 
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5^ 



9/i 



(76) 



Combining this result with Eqs. ([74|) and (|5J) gives the 
Gibbs Thomson relation in dimensionless units 



(5c(0, s) 
Ac 



= —doK{s) — l3v. 



where do is the capillary length given by, 



do = 



(Ac)2(aM/ac)e,' 

Cc + cTi/) is the total surface tension given by, 



du 



K 



if. 



(77) 



(78) 



(79) 



\ du J \ du 

and P is the coefficient of kinetic undercooling given by 
1 r cr^ 



f3 



(Ac)2(9^/9c) 



(80) 



Equation (|7^) provides a boundary condition at the in- 
terface for the diffusion equation of Eq. (^^ which can 
be written as 



d5c 
dt 



where 



D =T ^ 
^- ' dc 



This must be solved in conjunction with Eq. 
may be written, 



Acw(s) 



D, 



dSc 
du 



D, 



dSc 
du 



(81) 

(82) 
3|) which 

(83) 



0+ 



To understand of the significance of each term that 
enters the above equations it is useful to consider some 
limiting cases. First consider the case in which the con- 
centration field is a constant slightly different from Ceq, 
the equilibrium value, c = Ceq + Sc. If Sfi is the chemical 
potential difference between the phases defined by the 
non-conserved field at m = +00 and —00, then Eq. ( |7^ ) 
reduces to the AUen-Cahn equation in a field. 



Sn Ac 



(84) 



From this point of view the kinetic undercooling can be 
thought of as the simply the relaxation of surface ten- 
sion in a non-conserved field. Thus the Gibbs Thomson 
equation is equivalent to the Allen Cahn equation in the 
appropriate limit. 

The other simplifying case is when the non-conserved 
field is a constant as in a pure liquid or solid phase. In 



this case, the sharp interface equations remain the same 
except the coefficients a and /3 become 



and 



P = -i£/Ac)y{dfi/dc) 



(85) 



(86) 



For the conserved case of Model B, /3 is always negative. 
This term takes into account the lag of the concentration 
field behind a moving front. When the interface is mov- 
ing, the interfacial profile cannot instantaneously relax 
to the correct equilibrium shape Cq". For the conserved 
field, this correction is roughly as important as dynamic 
relaxation in the bulk phase, as will be seen in the next 
section. 



VI. LINEAR ANALYSIS 

To illustrate the influence of the various terms that 
enter the sharp interface model it interesting to study 
the dynamics of fluctuations around an almost planar 
interface. To flx ideas it is useful to consider an interface 
separating two phases that is defined by the equation y = 
h{x, t) as shown in Fig. (||). In the calculations to follow 
it will be assumed that V/i(x,i) is a small parameter. 
This is an additional constraint not implicit in the sharp 
interface models. 



u = {y-h(x,t)) cos(0) 




y-h(x,t) 
v'i+(ah/ax)2 



h(x,L) 



Phase I 



FIG. 6. Interface in Cartesian coordinates 

To facilitate the analysis it is worth noting that the 
normal velocity and curvature can be written in terms of 
derivatives of h as follows: 



du 1 dh 



(87) 
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and 



1 

7^ dx^ ' 



where 7 = -y/l + {dh/dxY 



A. Non-conserved dynamics 

As discussed in the preceding section, the sharp inter- 
face equations reduce to the AUen-Cahn equation in a 
field when the conserved field is a constant. When a sin- 

t valued interface as described above is considered, Eq. 
) reduces to the Kardar, Parisi, Zhang (KPZ) equa- 
tion 1 21] in the absence of noise. Substituting Eqs. ( |87| ) 
and (88) into Eq. ( ^4[ ) and Unearizing in h gives 



(89) 



dh _ d^h \ /dh^ 

where h —^ h—Xt, v = T^K^ and A = —vb^i^cja^. The 
additive noise term r\{x, t) in the standard KPZ equation 
|^l| appears when stochastic noise is included in the fun- 
damental Langevin equations. 

As a specific example, consider the following free en- 
ergy: 



« ,2 
-2^ 



+ -ip'' + d5cijj + -{dcY 



(90) 



where Sc = c — Ceq and c is a constant. For this free 
energy, a planar interface is stationary when d — s- 0. This 
interface is defined by the equations: 



ip'o = ^eq tanh 



^J■eq = 0, 



w 



(91) 



where ^ = K^j2a and V'eq = yajb- Thus the miscibil- 
ity gap, surface tension, and deviations of the chemical 
potential are given by 



Ac = AV'--2-V'e9, 

w w 



2 K^i'eq 



and 



5fi = -^Sc = wSc. 
oc 

Thus the coefficient A is given as 

^^ijj 3 

A = Tt(,K^{dSc) — r^(d(5c)- 

(Jt/) a 



(92) 



(93) 



(94) 



(95) 



B. Conserved dynamics 

Now consider an almost planar interface separating two 
phases of different concentration with the same free en- 
ergy as occurs for example in spinodal decomposition. 
Since concentration is a conserved field Eqs. (|77|), (|8l]), 
and (^3|) must be solved simultaneously. For simplicity, a 
two sided model in which dpi/dc is the same on both sides 
of the interface will be considered. This implies that the 
parameters do, P and D are the same in both phases. It 
is straightforward to perform the calculations in the more 
general case but this does not introduce any new physics 
and is not very illuminating. In the limit (dh/dx)'^ <C 1 
it is convenient to seek solutions of Eq. (pVj of the form 



Sc{u,x) = (5c(0)e(''=^-«l"l+"*\ 



(96) 



where u k, y — h{x,t). For this perturbation it is easy 
to show that the dimensionless quantities (D = ujd^/Dc, 
k = kdo, q = qdo and the dimensionless kinetic coefficient 
/3 = (3Dc/do satisfy 



-2f 



1 - 2q{l 



(97) 



-10 



-15 




-4-2 2 

log® 

FIG. 7. Dispersion relation for a planar interface sepa- 
rating two conserved phases of equal free energy. The four 
lines plotted in this figure correspond from right to left 
-f3 = 0, 10, 100 and 1000. 

In the long wavelength limit (k ^ 0), k ^ q and ui « 
—2k^, as expected. It is also interesting to note that in 
the short wavelength limit (k — > 00), q 1/(2 + 2\(3\) 
and Lo ~ — fc^. The crossover from the diffusion limited 
Lu ~ — fc^ at short wavelength to the asymptotic long 
wavelength behavior uj ^ —k^ occurs at smaller values 
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of k as the kinetic coefficient (3 becomes more negative. 
This behavior is sketched in Fig. (0). 

This analysis shows that the term 1 + (3 gives rise to 
diffusive (i.e., uj ^ —k^) behavior and is associated with 
relaxation in the bulk and in the interface shape. For 
example if solutions of the form g^'^'^-iM are sought in- 
stead of Eq. ( p6| ) this term becomes simply f3. Thus the 
'1' represents diffusive relaxation in the bulk and the /3 
represents relaxation of the interface shape. 

A model commonly used to study spinodal decomposi- 
tion is known as the Cahn-Hilliard model or Model B in 
the Halperin and Hohenberg classification scheme. The 
free energy for this model can be written: 



and 



dc 



iVcl' 



— = rcV2— . 



(98) 



(99) 



dt " 5c 

For this model a stationary planar interface is given by 



Ceq tanh 



where ^ = Kc/2a and Ceq — \f t 
entering the sharp interface equations are then: 

do ^ \i. 



(100) 



The coefficients 



/3=-6, 



and 



(101) 



(102) 



(103) 



C. Non-conserved and conserved Dynamics 

Now consider the stability of a stable phase invading 
a super-saturated liquid phase at constant velocity. This 
is precisely the situation considered by Langer et al ||] 
in the absence of kinetic undercooling and without relax- 
ational kinetics in the bulk phases. It is easy to show 
that the only solution for a planar front moving at con- 
stant velocity which is consistent with the sharp interface 
model (i.e., Eqs. (IT^), (|i|), and (|8|)) is 



Ac 



exp(-t;?;7Dc) - 1 



liv y' >0 

y'<0 



(104) 



where v is the velocity of the front and y' = y — vt is a 
coordinate in the co-moving reference frame. 



The stability of this moving front can be determined by 
studying perturbations about the planar front solution of 
Eq. (104). We seek solutions of the form 



5c_ 



Sep 
Ac 



(5;exp(i/c ■ X + utt — qy') y' > 
Ssexp{ik ■ X + ujt + q'y') y' < 



(105) 



where q,q' > and the position of the perturbed front is 
at 



y' — h{x, t) = /ifcexp(jfc • x + ujt). 



(106) 



To linear order in hk and Si^s, it is straightforward to 
show that the dimensionless w and k are determined by 

q'=q-2/l 

_ _ 2{q-2/l)il/l-q{q-l/l)) 



1-2(9-1/0(1-/3) 



^q{q-2/l)-LU 



(107) 



where I = 2Dc/doV and (3 = fiDc/do- The dispersion 
relation for u](k) is plotted in Fig. (H) for I = 100 and 
several values of [3. Note that, in contrast to Model B 
where the kinetic coefficient has a definite sign /? < 0, in 
Model C it c an h ave either sign. Note also that, when 
oo, Eq. (107) reduces to the result of Eq. (97) for 
conserved Model B dynamics, as it should. 




0.15 



FIG. 8. The linear dispersion relation for the MuUins Sek- 
erka instability including the kinetic coefficient (3, for J = 100 

A simplified version of the general system discussed 
in this paper has been extensively used to study single 
phase solidification phenomena [ p^ . The free energy T 
can be written as 



(108) 



where $ — c + h^ip) with 
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(109) 



With this form of /(V'), the interface width is W. The 
dynamical evolution is governed by Langevin equations 
for the conserved field c and the non conserved "0 which, 
in the noiseless limit, are 



'dt 

dc 



Dr 



1 5T 

T Sip 



dt bX Sc 
A stationary planar interface is given by 



(110) 



tanh 



WV2 



(111) 



For an interface with curvature k propagating with veloc- 
ity u, it is tedious but straightforward to use the Gibbs- 
Thomson relation of Eq. (|77|) to find 



/W 



4^/2 



.4\/2^ V WaJ ~ V 840 / v"d / 



(4^/2)(l 



(112) 



The term inside the square brackets of Eq. (112) is pro- 
portional to the kinetic coefficient P and contains the 
sum of the positive Model A contribution and the nega- 
tive Model B part. In principle the kinetic coefficient /3 
can be of either sign in solidification processes while it 
must be negative for any process described by Model B, 
such as phase separation in binary alloys. 
It is also trivial to show that 



dc „, 
at 



/ dc 


dc 




\du 


+ du 


J 



(113) 



These results are identical to those found numerically by 
others Pi. 



VII. SUMMARY 

The use of continuum phase-field models to describe 
phenomena involving the motion of well-defined sharp 
interfaces is discussed. The phase-field models involve 
interfaces which are diffuse on a length scale of ^. Con- 
sidering a general class of phase-field models, it is shown 
how equations describing the sharp-interface limit are ob- 
tained when <C 1 and ^v/D <C 1. It is also shown that 
the AUen-Cahn equation is a special case of the Gibbs- 
Thomson relation. 

In particular, it should be emphasized that these cal- 
culations are independent of the specific form of the 



free energy functional, provided describes well defined 
phases. Furthermore the calculations are universal: a 
large class of free energies give rise to sharp interface 
equations which differ only in the values of parameters 
but are of the same functional form. To realize this goal, 
it is essential that the "sharp-interface limit" involves an 
interface of finite width, ^. Expansions involving a zero 
width interface require a delicate and unphysical tun- 
ing of parameters in the free energy for thermodynamic 
consistency. In this work, the small parameters and 
vS^I D vanish when the curvature k and the inverse diffu- 
sion length V / D go to zero for a finite interface thickness 
^. Thus, delicate tuning is not required for thermody- 
namic consistency in our approach which is based on the 
fundamental principles of statistical mechanics. 

This work opens the way to construct physically con- 
sistent sharp interface descriptions of more complicated 
multiple phase systems such as a solid in contact with 
a fiuid which can support flows. This will involve mode 
coupling terms in the dynamical equations. Once such 
Langevin equations are constructed, there should be no 
conceptual difficulty in deriving the corresponding inter- 
face equations. 
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APPENDIX A: CURVILINEAR COORDINATES 

The curvilinear coordinates {u, s) used in the text are 
related to the Cartesian coordinates in the following man- 
ner 



r = R{s) + un(s). 



(Al) 



where R is the position of the interface and n(s) is the 
normal vector (see Fig. ^ . The metric tensor ga/s of the 
transformation from Cartesian to curvilinear coordinates 



IS 



1 

^~ I (1 + mk)2 



where k ~ dO/ds is the curvature with 9 the angle be- 
tween the X-axis and the tangent to the curve. The 
Laplacian in (m, s) is then obtained in the usual manner 



E 

Q,/3 



1 



d 



dx° 



'\9 



d 
d^ 
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K d 



1 



52 



&^ 

dv?' 1 + UK du (1 + uk)^ ds^ 

UK, d 



(1 + ukY ds 
inverse of the matrix g and = dn/ds. 



(A2) 



where a; = u, x'^ = s, g'^^ are the components of the 
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FIG. 9. Curvilinear Coordinates 

In the inner region, the fields vary much more rapidly 
in the u direction than the s direction. The coordi- 
nates {u, s) are rescaled in dimensionless units as (u, s) = 
(u/£,,es/^). The dimensionless curvature R = ^K/e and 
= In terms of these dimensionless vari- 

ables, the Laplacian becomes 
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where 
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= Rd/du 


(A5) 




= aVas^ - R^ud/du 


(A6) 




= -2uRd^/ds^ + R^u^d/du - uRgd/ds. 


(A7) 



APPENDIX B: GREEN'S FUNCTIONS 

It will also be useful to develop an expansion for the 
inverse of the Laplacian or Green's function. The Green's 
function of interest is defined by 



V|G(r,f') =(5(f-f'). 



(Bl) 



An expansion of the Green's function can be obtained in a 
straightforward manner. Let G{r, r ') = Gq{u, s; u', s') + 
eGi{u, s; u', s') + ■ ■ ■ where 

CoGo = (B2) 
CoGi + CiGo = S{u ~ u')S{s - s') (B3) 
C0G2 + CiGi + £2Go = (B4) 

and so on. The solution for Go is Go = 0, so that the 
lowest order solution for G is Gi , which satisfies the equa- 
tion 
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Gi {u, s]u',s')^ S{u - u')S{s - s') (B5) 



du^ 

which has the solution 



G^ {u, s; u' , s') 
G^(u, s; u' , s') = 



I -1^ _ { +u5{s - s') -C, <u' <u <Q 



+u'5{s - s') -( < u < u' < 0, 
■u'S{s-s') 0<u'<u<+C 



-uS{s - s') < u < u' < -fC- 
The surface terms of Eq. (|^) can be expanded as 



± / du^ i ds' 
lo ou Jb 

_(_ dSfi{u', s') 



G 



du' 



f^'^du^ i ds'yye^+^^-' 

Jo du Jb ±:1,±1 



/o JB 



n—1 7n—0 



dGt{u,s;u',s') 

Sf^n{u,s) — 

- G^{u,s;u,s) 



du' 



(B6) 



where 
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(B7) 



where the subscript B indicates that the integrands are 
evaluated on the boundary at u' = 0''' and at u' — 
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The 0{e) surface contribution becomes 



Si = 5r + 5+ 

du \ 

du \ du' 



du' 



-C 



du' 
dSn^ 



du' 



r+C 



duKi-o^c^m 



duK{+0-c],-{u)]. 



(B8) 



To evaluate the bulk contribution B, Eq. (|39| ) is ex- 
panded in powers of e. Eq. (|39|) reads 



c Jo 



dr"G±(f,f') 



9c'"(f',i) 



v± 



dt 



(B9) 



We note that dc"'{r',t)/dt = vdc"'/du' = ©(e) since 
the normal interface velocity v = evi + • • •, c™ = Cq" + 



e<5c'i° + • • • and G± = eGf 



Naive power counting 



seems to imply that — 0{e ), but changing variables 



f'^u', s' yields 



rrJo j^^i + u'eR{s')^^^^ 
= tBf + e^^^ + • • • 



(BIO) 



where 

Bf = I du I du'ds'viGf 
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(9u 9u' 5u du' 

For the Green's function introduced above, Bi becomes 



Bi = B^ +Bt = 



duK{u)-c°-\u)f 



(B12) 



APPENDIX C: TWO-SIDED MOBILITY 

In this appendix the sharp interface equations are out- 
lined for a mobility that takes on a constant value in each 



phase. There is the question of consistency of such a the- 
ory in the presence of stochastic noises in the underlying 
Langevin equations which we do not attempt to answer. 
We consider the system in the unphysical limit of zero 
noise. For convenience, the mobility in the phase in the 
region u < (m > 0) is denoted iXt)- 

When the mobility Fc is different in the two phases, 
the equation of motion for the concentration c becomes 



V Fc • V6^i 



dc ^/ 
— V Fr • V — 

dt \ ^ 6c 

= r^V^Sfi + (VFc) • VSfi. 
The procedure outlined in the main text gives 



(CI) 



v± 



,G±(r>-) dc' 



J dS' ■ (VV'G± - G±V'V 



-f / df'G^[r,r')\/5^l' 



F' 



(C2) 



The last term in Eq. (C2) is 0{e') and can be neglected. 
I t is st raightforward to repeat the calculations of Section 
( IV A ) for the velocity of the interface and for the chemi- 
cal potential at the interface. All results remain the same 
except, the diffusion constant Dc has an obvious depen- 
dence on FJ and 



Ac 



du 



(u) 



4" {u)] 

f ^[cjr(«)-cr(^)]' 

J — oo ^ c 



(C3) 



and the interface position m = is determined by the 
solvability condition 

r ^ [cr (u) - ci," (u)] = f ^ [4" (u) - cr in)] 

Jo i c J-oo t c 

(C4) 
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In two dimensions the tangential component of the time 
derivative can be written vrd/ds, where vt accounts for 
flow of concentration (or ip) along the interface, vt must 
be 0{e) since it is zero when e = 0. Furthermore in the 
inner region s is scaled with e so vrd/ds ~ 0{€^). When 
this term operates on c or ip is becomes 0{e^), since both 
c and i/) are independent of s. The time derivative in the 
moving frame is also at least 0{e^). This term 

accounts for the relaxation of fluctuations around the 
steady state profile not due to motion of the interface. For 
the concentration field fluctuations relax as t ~ (length)^ 
which implies ^ 0{e^). As with the tangential 

component this term becomes 0{e^) when operating on 
c. For fluctuations relax exponentially fast and thus 
this term can be neglected. 
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